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Abstract 

Generation of an acoustic wave by two pump sound waves is studied in a three-phase marine 
sediment that consists of a solid frame and the pore water with air bubbles in it. To avoid shock- 
wave formation the interaction is considered in the frequency range where there is a significant 
amount of sound velocity dispersion. Nonlinear equations are obtained to describe the interaction of 
acoustic waves in the presence of air bubbles. An expression for the amplitude of the generated wave 
is obtained and numerical analysis of its dependence on distance and on the resonance frequency 
of bubbles is performed. 
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I. INTRODUCTION 


Marine sediments as porous media consisting of a solid frame and the pore water are 
known to exhibit strong elastic nonlinearity M- In the present paper nonlinear sound 
wave generation by two pump acoustic waves is studied in a three-phase marine sediment, 
that is in a sediment that contains air bubbles in its liquid fraction. Nonlinear oscillations of 
gas bubbles in fluids lead to considerable enhancement of the nonlinearity of a homogeneous 
medium [5H8] as well as of the nonlinear properties of marine sediments that contain bubbles 
in the pore water 0 El] . In Ref. nu nonlinear interaction of three sound waves was studied 
in a two-phase marine sediment without air bubbles in the liquid phase. In the present paper 
the sound wave generation in a three-phase sediment is considered, as in Ref. HU, in the 
frequency range where there is a significant amount of wave velocity dispersion [IMS]. Due 
to velocity dispersion no shock-wave formation would occur which allows observing more 
explicitly interaction of three waves propagating in three different directions. 

II. THEORY 

Our task is to find the amplitude of an acoustic wave (cu 3 , k 3 ) generated by two pump 
acoustic waves (ca^ki) and (tJ 2 ,k 2 ), the nonlinearity being due to oscillating bubbles. The 
frequencies and wave vectors satisfy the energy and momentum conservation laws, 


uq + UJ2 — ca 3 , k x + k 2 — k 3 . 


The dependence of the wave (cu 3 , k 3 ) amplitude on the resonance bubble frequency and 
on distance will be analyzed along with the possibility of experimental observation of the 
process. 

To solve the problem we are based on the well-known Biot model in the form of the 
equations of continuity and momentum conservation for the liquid and solid fractions [ T6] . 



dv dp 

dt dx ’ 



( 1 ) 


In these equations we can limit ourselves for simplicity with one-dimensional case since 
the angles between the three waves are small enough. Indeed, numerical estimates HU 
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show that for the experimental situation which is used here and which is described in Refs, 
pans! the angles between the propagation directions of the waves (oq, k^), (cu 2 , k 2 ) and the 
propagation direction (let it be the x-axis) of the wave (u; 3 ,k 3 ) are approximately 9i ~ 25° 
and 0 2 ~ 18°. 

In Eqs. Q nonlinear hydrodynamic terms are omitted since as it was noted in Introduc¬ 
tion nonlinear acoustic processes are governed mainly by nonlinear oscillations of bubbles 
contained in the pore water. In Eqs. 0 p m is the mean density of the mixture of water 
and air bubbles; p s is the solid phase density, (the subscript ” 0 ” refers to the equilibrium 
values); v and u are the velocities of the liquid and solid phases; p is the pressure in the 
water; m is the porosity; a xx denotes the effective stress in a porous medium (see 1 (TI). 

°xx = - [k + (4/3 )p] — + p , 

POs rv s 

k and p are the bulk and shear moduli of the frame of the porous medium. 

Further calculations performed by one of the authors in Ref. m lead to the equations 


d 2 p m d 2 p 
dt 2 


v 


d 2 p s 


mn 


d 2 V 


d 2 p s 


(1 — m) — 


p of G dx 2 p 0s G dt 2 
k + (4/3)/i d 2 p s 


— v- 


G dt 2 
d 2 p 


= 0 


( 2 ) 


dt 2 p 0s dx 2 dx 2 

Here n is the bubble concentration and V is the bubble volume. Eqs. 0 should be supple 
rnented with the equation for the individual bubble motion [7], 


V + lu 2 0 V + fV - aV 2 - 2VV + V 2 ) = ep, 


(3) 


where ujo is the resonance bubble frequency. The coefficients in Eq. 0 are expressed through 
the equilibrium bubble volume Vo, its radius Rq and the adiabatic index 7 , 


a = Uq(1 + 7)/2Vb, fi = I/6V0, e = 4vrR 0 /po/, / = 


5 is the dimensionless absorption coefficient of bubble oscillations. 

As was mentioned in Introduction in order to observe three-wave interaction in a more 
explicit way we should use the frequency range with a noticeable velocity dispersion. It 
is preferable to choose the interval of the maximum velocity dispersion to ensure that the 
angles between the waves are not too small, otherwise they would fall in the dissipation 
spreading of the waves. As in Ref. m, the data from H 2 J listed in Fig. 2 of [H] will 
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be used. Basing on these data we choose the frequencies equal to uq = 2tt ■ 2 • 10 3 s _1 , 
cu 2 = 2vr • 3 • 10 3 s _1 , for which the sum frequency equals 2tt ■ 5 • 10 3 s _1 . We shall seek for 
the solution to Eqs. ([2]) , (|3]) in the form 

Pi = l/2[Pie i{kiX - Uit) + c.c.], Vi = 1/2 [y ie i(fci *” w<t) + c.c.], i = 1,2,3. 


The only nonlinear term in Eqs. ([2]) that ensures the nonlinear process is the last term 
of the first equation. In fact this term is the sum of the linear V 1 and the nonlinear part 
V n . The nonlinear part V n responsible for the generation of the wave p 3 by the waves p \, 
P 2 can be obtained from Eq. ^ in the form 

V n = _ [a - P(uj + u% + ^i^ 2 )]e 2 PiP 2 _ m 

(-^3 + w o - iuj 3 /)(w? - + iuj i f)(u% - UJo + iuj 2 f) ‘ 

In solving the equations ([2]) the so called approximation of a fixed held is used, that is 
the changes of the amplitudes Pi, P 2 because of nonlinearity are neglected, their changes 
being only due to linear bubble oscillations that cause dispersion in the medium. Using Q 
we obtain the solution to Eqs. ([2]) for the amplitudes P\, P 2 , Ps that are assumed to vary 
slowly at the distance of order of the wave length, 


dP 1 


= ib]P\, 


dP, 


— ib 2 P2i 


dx dx 

The following notations are used in Eqs. (J 5 ]) , 


dP 

-r - = ib 3 P 3 - iaPiP 2 - 
dx 


(5) 


bj = ( pfnVu j cK/2'yP 0 L ) (l - cuJ/wq - iSjUj/ufi) 1 , j = 1 , 2 ,3, 


where 


K — 1 — m — 


k + (4/3)/i 


p s c z 


L — K — —— 
m p s 


k + (4/3)// ^-r + 1 
p s c 2 


v = 1 — m — k/k s 


Po is the equilibrium pressure in the pore water, related to other equilibrium parameters of 
a bubble as Wq = 3yPo/ PqRq. The quantity a in the last of Eqs. ^ is the vertex of the 
nonlinear interaction, it equals 


_ p f nVu 3 cK [1 + 7 - (l/ 3 cUo)(^i + + OJ1UJ2)] 

47 2 P 0 2 L(1 - u\/ul - i5iul/ul){l - ul/ul - z<5 2 u;^/a;g)(l - cuf/cug - iS 3 ujl/u%) 

As can be seen from Eqs. (|5]) the imaginary parts of the quantities bj define dissipation of 

the acoustic waves. This dissipation is caused by energy losses due to bubble oscillations. It 

grows significantly when the resonance bubble frequency cuq approaches one of the frequencies 
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of the interacting waves. Numerical estimates show that in the frequency range where Uq 
is either much higher or much lower than the frequencies of the acoustic waves the wave 
attenuation caused by bubble oscillations is considerably less than acoustic dissipation of 
hydrodynamic nature. This dissipation does not enter Eqs. ([5]), but it is to be taken into 
account to estimate the real distance of the nonlinear interaction. The distance cannot 
exceed the dissipation lengths of the waves (aq, ki) and (u > 2 , k 2 ). The amplitude attenuation 
coefficients of these waves are equal correspondingly to aq ~ 0.8 ■ 10“ 3 cm -1 and to a 2 ~ 
3 • 10 _3 cm _1 (see Ref. |15j). This corresponds approximately to propagation distances 
~ 1250 cm and ~ 330 cm which means, that the interaction length does not exceed the 
distance that is the order of 300 cm. 

If at x — 0 the amplitude |P 3 | of the generated wave is at the fluctuation level, that is 
practically zero, the solution to Eqs. (|5]) can be represented in the form 

|P 3 | 2 = o^|Pi| 2 |P 2 | 2 P [1 + 7 - (l/3wo)O^i + 4 + (UKU,)] 2 /4 7 2 P 0 2 P, (6) 

with 

E(x) = e ~ 21mb3X + e ~ 21m ^ +b ^ x - 2e~ 21m ^ +b2+b3 > cos[Re(6 1 + b 2 - b 3 )x], 

D = (n 2 @ 2 @3 + (U2 ©i©3 + w 2 ©i© 2 + 2cnio; 2 03$i 2 — 2u}iU] 3 ®2&i 3 — 2a; 2 tt;30i < f ) 2 3. (7) 

In the formula (JT]) the following notations are used, 

©j = a 2 + jdj, Tjfc = a i a k + Pj/3 k , j, k — 1,2, 3, 

where 

= 1 - 4/4, 4 = s>4/4- 

We shall analyze the obtained amplitude of the generated wave, namely its dependence 
on distance at different resonance bubble frequencies and on the resonance bubble frequency 
at a given distance. To perform this analysis and for numerical estimates the following 
quantities and experimental data will be used PHHSII]EE], 

P 0 ~ 3 x 10 6 dyn/cm 2 , P 1 = P 2 « 10 5 dyn/cm 2 , 7 = 1.4, nV = 10 -5 , 

= 47t x 10 3 s _1 , cu 2 = 67T x 10 3 s _1 ; a; 3 = 107T x 10 3 s _1 , Sj ~ 4 x 10~ 2 , 

Pf = lg/crn 3 , p s = 2.65 g/cm 3 , m = 0.4, c ~ 1.7 x 10 5 cm/s, 
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k = 10 9 dyn/cm 2 , /i = 5x 10 8 dyn/cm 2 , k s = 3.6 x 10 11 dyn/cm 2 . 

Since the obtained formula Q is somewhat complicated computer simulation will be 
applied to investigate the behaviour of the sound-wave amplitude |P 3 | for various resonance 
bubble frequencies and various distances. 

1. The dependence of the ratio |P 3 |/|Pi| on distance for the resonance bubble frequency 
cuo that is less than the frequencies of the three waves, let it be cuo = 47 t x 10 2 s - 1 . The 
distance changes in the interval 100 cm < x < 300 cm. In this case the value of |P 3 | remains 
at any distance several orders of magnitude less than |Pi|. 

2. The resonance bubble frequency equals one of the frequencies of the three waves. The 
dependence of the ratio |P 3 |/|Pi| on distance in the interval 100 cm < x < 300cm is the 
following, 

a) u>o = oj\. The ratio |P 3 |/|Pi| is always rather small and is equal to ~ (3.6 — 3.5) x 10" 3 . 

b) cuo = CU 2 - The ratio |P 3 |/|Pi| becomes higher, of the order of « 3 x 10 ~ 2 and declines 
a little with distance because of dissipation due to bubble oscillations. 

c) lo 0 = l a 3 . The ratio |P 3 |/|Pi| is a little higher than in the two preceding cases, it is of 
the order of fa 4 x 10“ 2 , slightly decreasing with distance. 

The decrease of the ratio |P 3 |/|Pi| with distance in the last three cases is associated with 
the fact that the resonance bubble frequency in these cases always falls in resonance with 
one of the interacting waves which leads to enhancement of the energy losses due to bubble 
oscillations. On the other hand this evidences as well that in this frequency range nonlinear 
interaction is too weak to overcome dissipation. 

3. The resonance bubble frequency cuq is higher than the frequencies of the three waves, 
let it be equal to 107T x 10 4 s -1 . The distance changes in the interval 100 cm < x < 300 cm. In 
this case the amplitude |P 3 | grows slowly starting from the value ~ 0,17|Pi| up to the value 
~ 0.3|Pi|. This means that nonlinear oscillations of bubbles with high resonance frequencies 
ensures an effective generation of the sum frequency sound wave. 

In Fig. 1 the dependence of |P 3 |/|Pi| on the ratio cu 0 /cu 3 at the distance x = 300 cm 
is presented for the whole range of Uq starting from the frequency much lower than the 
frequencies of the interacting acoustic waves up to the frequencies that are much higher 
than the frequencies cui, 0 J 2 , and cu 3 . For cuo just between the frequencies uj\ and 0 J 2 there is 
a narrow peak that can be attributed to the fact that at this point loq is equally far from the 
resonances with Ui and cj 2 - For cjq > ui 3 there starts practically steady rise of the amplitude 
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FIG. 1: The dependence of |i- 3 1/1| on W 0 /W 3 at the distance x = 300cm (the interaction length). 


of the generated wave |P 3 | up to the value |P 3 | ~ 0.3|Pi 


III. SUMMARY 

The presented results show that the considered nonlinear acoustic process ’’feel” very 
slightly bubble oscillations with low resonance frequencies. Instead bubble oscillations with 
high resonance frequencies influence significantly nonlinear wave interaction. The generated 
wave amplitude in this case achieves a measurable value of the tenths of the pump amplitudes 
at a reasonable distance. The range where uo coincides with one of the frequencies of the 
interacting waves is not favorable for nonlinear acoustic interaction not only because of 
higher energy losses when bubble oscillations fall in resonance with the interacting acoustic 
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waves, but also because uq is not sufficiently high in that interval. 
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